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Long range polarization forces between ions and neutral atoms result in large elastic scattering 
cross sections, e.g., ~ 10 6 a.u. for Na-Na + or Na-Ca + at cold and ultracold temperatures. This sug- 
gests that a hybrid ion-neutral trap should offer a general means for significant sympathetic cooling 
of atomic or molecular ions. We present SIMION 7.0 simulation results concerning the advantages 
and limitations of sympathetic cooling within a hybrid trap apparatus consisting of a linear rf Paul 
trap concentric with a Na magneto-optical trap (MOT). This paper explores the impact of various 
heating mechanisms on the hybrid system and how parameters related to the MOT, Paul trap, 
number of ions, and ion species affect the efficiency of the sympathetic cooling. 



I. INTRODUCTION 

In 2003, W. Smith et ai. first proposed a hybrid ion- 
neutral trap consisting of a magneto-optical trap (MOT) 
concentric with and encompassed by a linear radio- 
frequency quadrupole (RFQ) Paul trap [HQ- Since then, 
other proposals have been made for the sympathetic cool- 
ing of molecular ions within similar hybrid ion-neutral 
traps 0]. Several experiments using hybrid traps have 
measured charge exchange cross sections for Yb-Yb + [H , 
Rb-Ca+ @, and Ca-Yb + @. Single Ba+, Yb+, and Rb+ 
ions have been sympathetically cooled within a hybrid 
Paul trap using a Rb Bose Einstein Condensate (BEC) 
to energies equivalent to sub-Kelvin temperatures (343 ■ 
Also, experimental evidence of sympathetic cooling of 
Rb + ions within a hybrid Rb MOT Paul trap was shown 
in Ref. [nj. 

Sympathetic cooling occurs when one gas is cooled by 
a colder gas via elastic, inelastic, and charge exchange 
scattering. Sympathetically cooling ions within a Paul 
trap is a more general technique than direct laser cooling 
[lll[T^|. since specific laser-excitable resonant transitions 
are not required of the trapped ion species. This tech- 
nique is useful for atomic species but is often the only 
option when trying to cool molecular ions [HI, [l4| . 

In ion-ion sympathetic cooling, one ion species is di- 
rectly laser cooled and then collisionally cools the other 
species. Due to the strong Coulomb interaction V oc 1/r, 
this is a highly effective method for cooling ions to cold 
and ultracold temperatures fl3l4T7l ] . The most general 
form of sympathetic cooling is neutral buffer gas cool- 
ing. A lower bound on the cooled ion's equilibrium tem- 
perature is set by the temperature of the neutral buffer 
gas. Additionally, the technique works best for ions 
(with mass mj) and neutral atoms (with mass m n ) whose 
masses meet the criterion mj/m n > 1, else the ion trap's 
inherent atom-ion rf heating mechanism can overwhelm 
the collisional cooling [18l - [22l | . 

In a hybrid trap, cooling by the MOT or BEC acts 
as a combination of the two previously mentioned tech- 
niques. The neutral species (MOT or BEC) is directly 



laser cooled, but it also acts as a highly localized cold 
or ultracold buffer gas. Due to the laser cooling and 
trapping of the neutral species, ions overlapped with a 
MOT or BEC could reach lower final temperatures than 
if they were overlapped with either a room tempera- 
ture or chilled buffer gases. Unlike buffer gas cooling, 
we show that a MOT can efficiently cool equally mas- 
sive ion-neutral species. Equal mass ion-neutral cooling 
has been observed experimentally within hybrid traps in 
Refs. UM- 

In contrast to ion-ion sympathetic cooling, 
it has been theorized that a hybrid trap should simulta- 
neously cool internal degrees of freedom as well as the 
translational motion of molecular ions [U, Q • 

R. Cote et al. have shown that the elastic scattering 
cross sections for both Na-Na + and Na-Ca + are large (~ 
10 a.u.) when compared to neutral-neutral or ion-noble 
gas (neutral buffer gas) cross sections in the relevant 
temperature regime (10 -3 to 10 3 K) @, El, H3] . This is 
due to the long-range polarization potential V oc — a/r 4 , 
where a is the dipole polarizability of the neutral species. 
These large elastic scattering cross sections suggest that 
the hybrid trap should offer significant sympathetic cool- 
ing. 

Using simion 7.0 software [HI, HI], we have simulated 
the sympathetic cooling of Ca + or Na + to energies equiv- 
alent to cold temperatures within a hybrid Na MOT and 
linear RFQ Paul trap. These custom simulations model 
experimental work currently underway in our laboratory. 
Simulations like these have proven vital to the under- 
standing of many Paul trap or hybrid trap experiments 
[l,[l6|, [20,113], where several papers have specifically used 
the SIMION software 0, [HI, SI [H, [H, SI] . Although the 
simulations presented in this paper model our actual hy- 
brid system, more general conclusions may still be drawn. 
We find that even in the case of modest MOT densi- 
ties (10 9 cm -3 ) and modest MOT temperatures (1 mK), 
single ions can be cooled to energies equivalent to cold 
and ultracold temperatures within a few seconds. These 
MOT conditions can also sympathetically cool more than 
one trapped ion although much higher MOT densities are 
needed to reach sub-Kelvin temperatures. 

This paper is organized as follows: In Sec. [Ill we review 
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the workings of our hybrid trap and the details of our 
simulation model. Next present the results of our simu- 
lations in Sec.lIIII We first discuss single ion sympathetic 
cooling ( Sec. IIII A[) and then multiple co-trapped ion 
sympathetic cooling, both in the ion cloud ( Sec. IIII B ip 
and ion crystal phase ( Sec. IIII B 2p . We conclude in 
Sec. IIVI Last, we include an appendix which discusses 
various inherent Paul trap heating mechanisms which 
limit the sympathetic cooling ability of hybrid traps. 



II. BACKGROUND 
A. Na MOT and linear rf Paul trap 

An illustration of the part of the hybrid trap inside 
our vacuum chamber is shown in Fig. Ufa). The Paul ion 
trap is comprised of eight end segments and four central 
rf segments (i.e., electrodes). Passing through the sides of 
the trap are six 589 nm MOT beams forming a standard 
Na MOT [3(||3l| concentric with the ion cloud that forms 
in the center of both traps. The MOT's magnetic field 
gradient is created outside of the vacuum chamber by two 
coils in an antihelmholtz configuration. 

Also shown in Fig. Ufa) is a 405 nm photoionization 
beam collinear with one of the MOT beams. The pho- 
toionization beam ionizes excited 3P 3 / 2 Na atoms within 
the MOT or the background Na gas. The ionization 
process is known as resonance-enhanced-multiphoton- 
ionization (REMPI) [32| and is one example of how ions 
can be loaded within the trap experimentally. 

For the remainder of this section we will focus on a brief 
review of Paul trap principles, terminology, and quanti- 
ties of interest within the context of our actual Paul trap. 
For a more detailed discussion of Paul trap physics see 

Refs. mm nisi. 

The four central Paul trap segments provide radial con- 
finement with an applied oscillating quadrupole driving 
field, giving the effect of a rotating saddle potential well 
or harmonic pseudopotential [35H38|. The driving field 
oscillates at angular frequency ft and has amplitude ±Kf 
(relative to electrical ground) on each diagonal pair of 
segments. The eight smaller end segments allow for axial 
confinement and are held at a dc potential V cn d during 
trapping. The unitless efficiency factor i] depends on the 
geometry of the end segments. 

All together, 



$(xi,i) ~ V T f cos(f^) ■ 
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is the total approximate (r <C r ) time-dependent electri- 
cal potential $ near the trap's center, with magnitude of 
each component of the position vector Xi [see Fig. HJb)], 
inter-electrode inscribed radius r , and rf segment length 
2zq. Near the trap's center, we can approximate the trap 




FIG. 1. (a) Illustration of the part of the hybrid trap system 
within the vacuum chamber. A Na MOT (orange) is over- 
lapped and concentric with an ion cloud (green) inside the 
segmented Paul trap with six 589 nm MOT beams (yellow) 
and one collinear 405 nm photoionization beam (blue). The 
photoionization beam ionizes excited 3P Na. (b) Axial view 
of Paul trap with Cartesian coordinate system 



electrodes' shape as hyperbolic and the vacuum chamber 
(at electrical ground) to be infinitely far away. 

The equation of motion for a single ion within the 
electrical potential described by Eq. (fTJ) is known as the 
Mathieu equation. For a single ion with charge e and 
mass mi, approximate solutions to the Mathieu equa- 
tion that are stable against ejection of the ion from the 
trap are possible for particular ranges of the so-called 
stability parameters dj and qi (ai < and < qi < 0.9) 
US HE H3, Hl| • The terms ai and qi are defined as 



ai = a 2 = — — 
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and q 3 = 0, <?i 
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The ion's motion within the harmonic pseudopotential 
can therefore be described as a superposition of slow sec- 
ular motion, with angular frequency 
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(such that ai and qi <C I) and micromotion at the driving 
field frequency f2 (whose amplitude increases as the ion 
moves farther away from the trap's nodal line) [34j . 

The total time averaged (denoted by ( )) kinetic energy 
(Ek) of the ion is defined as 



(E k ) = -mj (v?) = 



(4) 
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The secular motion determines the kinetic energy in the 
x 3 direction [34| • Therefore, the x 3 amplitude of the ion's 
motion 2 3 can also be expressed as a function of the ion's 
energy 
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assuming that mode approximately contains energy of 
where is the Boltzmann constant and T is the 
equivalent temperature associated with the ion's mean 
energy (E k ) = §fc b T 

It should be emphasized that the above discussion ap- 
plies to a single ion in a Paul trap under ideal vacuum 
conditions. As one introduces other ions or background 
gas collisions, the multi-body problem quickly becomes 
too difficult to solve analytically and the need for numer- 
ical simulations arises. 

When ion-neutral and ion-ion collisions occur, various 
ion heating mechanisms arise that are inherent to Paul 
traps. Within the scope of this paper, the two most im- 
portant mechanisms are atom-ion rf heating [l8l . [2lj and 
ion-ion rf heating [53, IH, El[ . Atom- ion rf heating 
occurs at certain instances of the driving field's phase 
when an ion's speed is instantaneously reduced by an 
ion- neutral collision resulting in a transfer of energy from 
the driven micromotion to the secular motion (2l| . Atom- 
ion rf heating should be differentiated from instantaneous 
collisional heating, which is when an ion-neutral collision 
simply increases the instantaneous speed of the ion (in- 
dependent of the ion's micromotion). Unlike a buffer gas 
which fills the entire volume of the Paul trap, a hybrid 
trap's neutral species occupies a finite region of the trap. 
As a result, certain approximations made in Ref. [l8[ 
do not apply and the restriction on equally massive ion- 
neutral sympathetic cooling imposed by atom-ion rf heat- 
ing within a buffer gas is not valid for a hybrid trap [l(| . 
If there is a high enough ion cloud density, ion- ion rf heat- 
ing results in an absorption of energy from the driving 
field due to chaotic motion within the ion cloud result- 
ing from the Coulomb interaction between the co-trapped 
ions [33| . A more detailed discussion of these mechanisms 
(and more) can be found in Appendix [Al 



B. SIMION simulation details 

SIMION 7.0 uses a fourth order Runge-Kutta method 
to numerically solve for the ion trajectories generated by 
the fields, produced by both the Paul trap's electrodes as 
well as the Coulomb repulsion between ions [III [26| . Our 
program simulates both the slow secular motion and the 
fast micromotion. We have created a custom electrode 
geometry which models the exact dimensions of our Paul 
trap, the vacuum chamber, and our ion detection appa- 
ratus (a biased mesh and Channeltron). Figure [5] shows 
SIMION 's cutaway 3D rendering of the electrode model 
as well as ion trajectories (blue) during ion extraction 
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FIG. 2. SIMIOn's isometric view of a 3D rendering of our 
segmented Paul trap after extracting ions for detection. The 
(blue) tracks show the ions' trajectory out of the center of the 
trap where 100% of the ions reach the Channeltron. 



and Channeltron detection. The axial ion extraction oc- 
curs (and can be simulated) when a gated dipole field 
is applied to the end segments. SIMION not only allows 
the user to build custom electrode geometries, but also 
contains a user programming interface that can be cus- 
tomized to control time-dependent fields, ion-neutral col- 
lision effects, and initial conditions [22| . 

Single ion simulations are always initialized at the cen- 
ter of the trap with an initial velocity azimuthal and polar 
angle of 45°. In multiple ion simulations, the ions' ini- 
tial spacial and velocity directions are isotropically dis- 
tributed. The energy of an ion at t — s is always set to 
the mean energy associated with the temperature of the 
neutral gas from which the ion is born. Typically this is 
from a 1000 K (0.1 eV) background gas. 

The program's time step At is continually adjusted 
such that the ion moves a specified number of grid 
units per time step (typically At ~ 10~ 3 — 10 _1 fj,s < 
rf period ~ 1 /is). The program simulates three envi- 
ronments: ideal vacuum conditions, a hot low-density 
neutral background gas, or both a background gas and 
high-density cold MOT. When running in either of the 
non-ideal vacuum environments, the probability of an 
ion-neutral collision is calculated within each time step 
according to 



Pa* = 1 
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where n is the density of the gas [19|, [22| . In Eq. ^ , the 
program uses cither the background gases' density or the 
MOT's density for n depending on the mode of operation 
and the instantaneous position of the ion (e.g., it uses the 
MOT density if the ion is inside a small sphere specified 
by tmot)- 



K S {E) = a s (E)v = a a (E) 




(7) 



4 



is the instantaneous rate coefficient associated with the 
atom-ion elastic (s = el) or non-radiative charge ex- 
change (s = ce) scattering cross sections as a function 
of the instantaneous collision energy E, relative velocity 
v, and reduced mass \x. The cross sections 
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with coefficients C c \ (4174 a.u. for Na-Na + and 5070 
a.u. for Na-Ca + ) and C ce (57 a.u. for Na-Na + ) were 
calculated using power-law fits from a quantal ab-initio 
treatment in Refs. 0, [24[. We find that the mean time 
between collisions depends on several parameters (most 
importantly neutral species density), but is typically in 
the range of ~ 10 2 — 10 4 /is. 

Using a random number generator and Eq. ^ , the pro- 
gram decides whether or not an instantaneous collision 
will occur during each time step [Til in the event 
a collision occurs, the neutral atom's initial speed and 
direction are chosen using a random number generator. 
The generated speeds adhere to the Boltzmann distribu- 
tion and the initial direction is isotropically distributed. 

The ion's final velocity during a charge exchange col- 
lision is determined by swapping the ion's current veloc- 
ity with the randomly generated velocity of the neutral 
atom. In an elastic collision, within the center of mass 
frame, the final velocity of the ion is calculated and forced 
to adhere to a pseudo-hard sphere differential scattering 
cross section The ion's azimuthal scattering angle is 
isotropically distributed and the polar angle 6 follows the 
distribution function p described by 
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The various randomly generated distribution functions 
[Boltzmann, isotropic, p(E, 6)] are created via a Monte 
Carlo select and reject method [421 ] . 

There is a precedent for usin g a hard sphere model in 
these types of simulations [lSj |43I445| . However, when 
similar systems have been analyzed with a full quantal 
treatment (e.g., Yb-Yb + ), the differential cross sections 
have not been found to be isotropic within the tempera- 
ture regime being considered (10 -3 to 10 3 K) (46[. This 
has also been observed experimentally [2l], [22j ■ The sys- 
tem must be in the nanoKelvin regime to exhibit pure s- 
wave scattering. A fully quantal treatment considers the 
higher order partial wave contributions, which generally 
results in a differential cross section that favors forward 
scattering. To improve upon the hard sphere approxima- 
tion, we still use a rectangular differential cross section, 
but it is only nonzero for angles less than 60°. This cross 
section is what we are calling the pseudo-hard-sphere 
differential cross section. When results were compared 
using a true isotropic hard-sphere differential scattering 
cross section, we found thermalization times to be slightly 
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FIG. 3. Plot of kinetic energy vs. time of trapped single ion in 
different program environments. Curve (a) is from a Na + ion 
within background Na gas only (black) Pback = 7 x 10~ 9 torr 
and T back = 1000 K (0.1 eV). Curve (b) is from a Na+ ion 
under ideal vacuum conditions (magenta). Curve (c) is from 
a Na + ion cooled by a MOT (orange) with n = 5 x 10 9 cm -3 
and Tmot = 1 mK (10 -7 eV) within a background gas 
(Pback = 10" 9 torr and T bac k = 1000 K (0.1 eV)). The curve 
is fit using Eq. {TO}. Curve (d) is from an initially cold heavy 
ion with mi/m„ ~ 0.26 (red) heated by a MOT under the 
same neutral gas conditions as (c). 



shorter, but final temperatures to be approximately un- 
changed. 

Unless otherwise specified, the energy values and root 
mean squared positions reported throughout this paper 
are a time average over ss 15 secular oscillations of the 
instantaneous kinetic energy of a single ion (or mean en- 
ergy of a group of ions) queried once per time step. 

Our simulation reproduces results consistent with the 
existing Paul trap literature. For example, under ideal 
vacuum conditions, with no excess micromotion, and us- 
ing optimal stability parameter settings, a single ion has 
no heating mechanism. Therefore, when a single Na + 
ion was simulated under ideal vacuum conditions, the nu- 
meric precision was increased until the ions mean energy 
remained constant [Fig. [3] curve (b)]. 

When interacting with only a hot low-pressure neutral 
background gas the ion (initially at the mean energy asso- 
ciated with the background gas's temperature) heats up 
due to atom-ion rf heating and instantaneous collisional 
heating. Additionally, a single ion has far fewer colli- 
sions with the background gas as compared to its interac- 
tion with both a background gas and a cold high density 
MOT. For example, 21 collision events with a background 
Na gas can be clearly seen as discontinuities in curve (a) . 
In contrast, curve (c) shows sympathetic cooling (from 
~ 0.1 eV to ~ 10~ 6 eV) after 256 elastic scattering colli- 
sions and 48 charge exchange collisions with atoms from 



5 




0.0 0.5 1.0 1.5 2.0 

Time (s) 



FIG. 4. Plot of axial position of one Ca + ion relative to 
the trap's center vs. time when overlapped with a MOT 
n = 6 x 10 9 cm" 3 , Tmot = 1 mK (10" 7 eV), and back- 
ground gas Pback = 10" 9 torr with T bac k = 1000 K (0.1 eV). 
As the ion is cooled the axial amplitude decreases. 

the modestly dense (n = 5x 10 9 cm -3 ) and cold MOT 
[Tmot = 1 mK (1CT 7 eV)]. We see that the hybrid trap 
can yield sympathetic cooling despite the high atom-ion 
rf heating associated with mi/m n « 1. If the mass ratio 
becomes m / / m n > 1 , atom ion rf heating is reduced and 
greater cooling can be achieved, as depicted in Fig. SI 
which shows a single Ca + ion cooled under similar MOT 
conditions. As the Ca + ion is cooled the axial oscillations 
approach zero amplitude in accordance with Eq. ([5]). If 
the mass ratio is mi/m n < 1, atom- ion rf heating colli- 
sions with the cold MOT can actually heat an initially 
cold single ion [seen in Fig. [3] curve (d)] as predicted in 
Ref. [l||. 

Since the energy dependence in K e \ is weak we can 
approximate the net heating and cooling rates to be a 
constant k. Therefore the time dependence of the ion's 
energy can be approximated as 

E(t) ~ -Eflnal + (-^initial ~ ^final)e _f "* , (10) 

which our simulated ion's energy evolution follows in 
Fig. [3] curve (c). 

III. RESULTS 

A. Single ion 

The settings for all simulations, unless otherwise spec- 
ified, are the following: the Na MOT is concentric with 
the Paul trap (where there is zero micromotion ampli- 
tude), Tmot = 1 mK (10" 7 eV), P bac k = 1 x 10" 9 torr, 
Tback = 1000 K (0.1 eV), V cnd = 35 V, V rf = 40 V for 
Na+ (or 70 V for Ca + ), fl = (2tt)708 kHz, r] « 0.14, 
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FIG. 5. Plot of kinetic energy vs. time for a single ion 
showing the effect of MOT density on sympathetic cooling 
[Tmot = 1 mK (10~ 7 eV) and tmot = 1 mm for all curves]. 
Curve (a) is from a Ca + ion cooled by a MOT with 
n = 1 x 10 9 cm -3 (blue). Curve (b) is from a Na + ion cooled 
by a MOT with n = 5x 10 9 cm -3 (orange). Curve (c) is 
from a Ca+ ion cooled by a MOT with n = 2.5 x 10 10 cm" 3 
(navy). Curve (d) is from a Na + ion cooled by a MOT with 
n = 2.5 x 10 10 cm" 3 (red). Higher MOT density results in 
lower final energy and faster thermalization. 



ro = 9.5 mm, and zq — 24 mm. Therefore, q\ rs 0.4 and 
03 ~ °' 2 m amu - These values were chosen to match closely 
with our optimal experimental settings and actual trap 
geometry. 

We have found that the cooling rate and final temper- 
ature of the ions depend on several parameters with the 
MOT density being the most critical, due to the exponen- 
tial n dependence in Eq. ([6]) . Collisions with background 
gas atoms at pressures below 10" 8 torr (easily obtain- 
able experimentally) had a negligible effect on sympa- 
thetic cooling. We observe only one or two background 
gas collisions out of hundreds or thousands of MOT atom 
collisions at these densities. 

As the MOT density increases, the thermalization time 
and equilibrium energy decrease. For example, a Ca + 
ion overlapped with a MOT of density 1 x 10 9 cm" 3 as 
shown in Fig. [5] [curve (a)] does not thermalize until 
~ 5 s and has a final energy of ~ 10" 6 eV, while Ca + 
cooled by a MOT with density 2.5 x 10 10 cm" 3 equili- 
brates at ~ 10~ 7 eV in ~ 0.75 s [curve (c)]. A single Na + 
ion shows the same trend, as can be seen by comparing 
curves (b) and (d) in Fig. [5J 

The jagged appearance of the curves in Fig. [5] can be 
attributed to the competing effects of instantaneous col- 
lisional heating, atom-ion rf heating, and instantaneous 
collisional cooling. The simulations show that the domi- 
nant heating mechanism is atom-ion rf heating. For ex- 
ample, only 7% of all elastic scattering collisions within 
the simulation associated with Fig.[S]curve (b) resulted in 



6 




FIG. 6. Plot of kinetic energy vs. time showing the effect 
of ion species on sympathetic cooling. Curves (a) Na + (or- 
ange) and (b) Ca + (blue) were cooled with a MOT having 
n~5x 10 11 cm -3 . Ca + cools to a lower final energy and ex- 
hibits fewer fluctuations than Na + , due to the reduced atom- 
ion rf heating. 



FIG. 7. Plot of kinetic energy vs. time for (a)-(b) single Na + 
(orange) and (c)-(d) single Ca + (blue) showing the effect of 
qi on sympathetic cooling and its dependence on ion species, 
(a) and (c) have the Paul trap stability parameter qi ~ 0.4, 
while (b) and (d) are at an increased V r { resulting in Paul trap 
stability parameter qi ~ 0.75. The larger atom-ion rf heating 
associated with qi ~ 0.75 overwhelms the MOT cooling. 



an instantaneous speed increase, i.e., instantaneous colli- 
sional heating. The infrequent number of instantaneous 
collisional heating events is likely due to the difference 
between the ion's equilibrium energy and the mean neu- 
tral atom energy associated with the MOT's tempera- 
ture. Further support for this explanation comes from 
the fact that if the ion's equilibration time and equi- 
librium energy is lower, the percentage of instantaneous 
collisional heating events increases [e.g., 20% within the 
simulation associated with Fig. [5] curve (a)]. Addition- 
ally, due to the reduced atom-ion rf heating associated 
with a larger mass ratio mi /m n we observe the following: 
smoother energy equilibration, lower equilibrium energy 
(e.g., 3 x 10~ 7 cV for Ca + compared to 3 x 10~ 6 eV for 
Na + in Fig. [5]), and smaller standard deviation energy 
fluctuation at equilibrium (e.g., 3 x 10~ 7 eV for Ca + 
compared to 2 x 10~ 5 eV for Na + in Fig. 

The atom-ion rf heating increases with increasing V T f 
due to the micromotion's dependence on the stability pa- 
rameter qi. By varying V T [ we found that absolute V r { 
values were not a good metric for atom-ion rf heating 
rate comparison between different ion species, but q\ was 
[38] . Figure [7] shows MOT sympathetic cooling of Na + 
in curves (a) (qi ~ 0.4) and (b) (gi ~ 0.75), while Ca + is 
shown in curves (c) (qi ~ 0.4) and (d) (q\ ~ 0.75). Un- 
der ideal vacuum conditions increasing q\ only caused a 
small increase in the heating rate, likely due to operating 
the trap close to the upper edge of the single ion stability 
boundary. Hence, the difference between the left plot and 
the right plot for a given ion species is almost entirely due 
to atom-ion rf heating. The sympathetic cooling cannot 
combat the heating from high rf amplitudes. Therefore, 
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FIG. 8. Plot of final thermalized energy of a single Ca + ion 
(blue) circles and Na + (orange) triangles vs. the MOT's tem- 
perature. As the MOT temperature is lowered the final mean 
energy of the ion decreases, but not indefinitely, since the 
atom-ion rf heating rate increases as well. 



it is necessary to use low q\ values (qi < 0.4), provided 
the trap depth is not lowered below the initial energy of 
the ion. Single ion experiments within hybrid ion-BEC 
traps have drawn similar conclusions 7\. For multiple 
ion cooling, using low q\ values will also eliminate any 
instability heating (l7| . 

Initially the ion's equilibrium energy decreases with 



7 



decreasing MOT temperature, but not indefinitely (see 
Fig. In fact, at low enough MOT temperatures the ion 
actually begins to equilibrate at higher energies. Again, 
at a given MOT temperature Ca + is cooled to a lower 
final energy than Na + . We believe this effect is caused 
by atom-ion rf heating. At lower MOT temperatures, the 
approximation made in [2l| (see Appendix that the 
neutral atoms have exactly zero velocity in the laboratory 
frame (causing collisions to result in large instantaneous 
speed decreases resulting in greater atom-ion rf heating) 
becomes increasingly more valid. Hence, at lower MOT 
temperatures the atom-ion rf heating rate increases, re- 
sulting in a higher equilibrium energy. In addition, we 
found that lowering the MOT's temperature had little 
effect on decreasing the thermalization time of the ion. 

The initial ion energy was varied as high as 0.7 eV 
(~ 3000 K), which resulted in little to no difference in 
final energy and thermalization time. R. DeVoe found 
similar results for buffer gas cooling of a single ion [2(| ■ 

We found that with a fixed number of Na atoms in 
the MOT (N = 5 x 10 7 ), a smaller MOT radius cooled 
faster and lower than a large MOT radius, i.e., increased 
MOT density is favorable despite decreased initial over- 
lap between the MOT cloud and the single ion trajectory 
volume (i.e., the volume occupied by the ion's 3D orbit). 
The ion's initial secular axial amplitude was ~ 2 mm (al- 
ways larger than the MOT radii tested, e.g., 0.25 mm- 
1.5 mm). This counter-intuitive result can be explained 
by the fact that with a higher density in the exponent 
of Eq. ([6]), there is a higher collision rate which is ap- 
parently more important than the reduced percentage of 
time spent initially overlapped with the MOT. 

To increase overlap without changing MOT character- 
istics, we compressed the initial ion trajectory volume 
by increasing the end segment voltage V cn d in Eq. (JTJ 
(although the initial amplitude was still larger than the 
radius of the MOT). We found that this offered little im- 
provement in thermalization time and final energy. Over- 
lap is improved automatically as collisions with the MOT 
cool the ion and decrease the ion's oscillation amplitude 
[in accordance with Eq. ((SJ) and seen in Fig. 0] . The fact 
that the ion's final energy is insensitive to MOT over- 
lap is consistent with the lack of sensitivity to the initial 
ion energy, given the connection between ion energy and 
secular oscillation amplitude described by Eq. (|SJ). 

For ions with no laser-excitable transitions, such as 
Na + , overlapping the MOT with the center of the ion 
trajectory volume becomes experimentally challenging, 
because there is no fluorescence to visually confirm ion- 
neutral concentricity. Therefore, we simulated the cool- 
ing of a single Na + ion by a MOT displaced axially off 
center, (see Fig. |SJ The MOT's ability to sympatheti- 
cally cool is dramatically reduced if it is not concentric 
with the ion cloud. If the ion tends to have collisions 
with the MOT at its secular oscillation turning point 
(when the ion is at its peak micromotion amplitude), a 
greater percentage of ion-neutral collisions will result in 
atom-ion rf heating. Also, as shown in Fig. SI as the 
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FIG. 9. Plot of kinetic energy of a single Na + ion vs. time 
showing the effect of MOT-Paul Trap concentricity on sym- 
pathetic cooling. Curves are from a MOT (pmot = 1 mm 
and n = 5.7 x 10 9 cm -3 ) (a) located 2 mm off center axially 
(green), (b) 1.5 mm off center axially (black), (c) 1 mm off 
center axially (purple), and (d) on center (orange). The ion's 
equilibrium energy is sensitive to reduced MOT concentricity 
greater than one MOT radius. 

ion cools its axial amplitude will decrease. Therefore, 
the ion's final energy is now limited (at best) to the en- 
ergy equivalent to the secular oscillation amplitude that 
equals the distance between the edge of the MOT and 
the center of the ion trajectory volume. This minimum 
ion-MOT concentricity amplitude is why the the final 
ion energy is not significantly affected until the MOT 
displacement is greater than one MOT radius, as seen 
in Fig. [9] A secondary consequence of the offset is a re- 
duction in overlap resulting in a smaller effective collision 
rate. Reference Q, using an ion-BEC hybrid trap, exper- 
imentally demonstrates that the effective collision rate is 
rather sensitive to ion-neutral concentricity which is in 
qualitative agreement with our findings. 

B. Multiple ions 

1. Ion cloud 

When simulating multiple ions, additional complexities 
were incorporated into the program. A single ion can be 
initialized at the center of a trap, but multiple ions must 
be distributed throughout space in an ion cloud. We ini- 
tialized the ions isotropically within a sphere concentric 
with the center of the Paul trap. In doing so, we found 
that the initial time averaged energy of the ion cloud 
(i.e., (Ek) of the ion cloud after at least one secular os- 
cillation) was highly sensitive to the initial size of that 
cloud and not the velocity of the ions within the cloud at 
t = s (see Fig. 10). The total energy of an ion at the 
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FIG. 10. Plot ol initial average kinetic energy (after two sec- 
ular oscillations) vs. (r rms ) 2 (at t = s) of a Na + ion cloud. 
Triangles (magenta) are for Ni on = 10 and circles (black) are 
for Ni on = 200. The initial kinetic energy of the ion cloud 
after two secular oscillations is almost entirely dependent on 
the (r rms ) at t = s and not the ions' 10 -7 eV (1 mK) ki- 
netic energy at t = s used for each data point. The energy 
is only approximately quadratically dependent since Eq. Q 
is cylindrically symmetric and not spherically symmetric. 



moment it is born is primarily determined by its large 
potential energy derived from its position relative to the 
trap's center, not its smaller kinetic energy derived from 
its pre-ionization neutral atom velocity. Therefore, the 
size of either the MOT or ionization beam (whichever is 
smallest) is what primarily determines the time averaged 
energy of an ion cloud after one secular oscillation. We 
concluded that ions created directly from a MOT inside a 
Paul trap will not have cold initial time averaged transla- 
tional energy, despite the fact that the neutral ensemble 
from which they are born is cold. 

The dependence on initial (r rms ) is also true for ions 
born from the non-localized background gas. Since ions 
born from the background gas may be born farther from 
the nodal line and have much greater initial velocities, 
ions born from the MOT will still be initially colder (but 
never cold or ultracold). 

Another level of complexity that needed to be consid- 
ered, once we allowed for multiple ion trapping, was the 
production of additional ions born directly from the Na 
MOT during the cooling period. Molecular ions NaJ 
created via photoassociative ionization and subsequently 
fast (~ 0.5 eV) atomic Na + created via photodissociation 
can be produced by the Na MOT's trapping lasers (47i — 
l50l |. Although Na is the only alkali that undergoes this 
photoassociative ionization from its own MOT beams, 
alkaline earth MOTs (e.g., Ca, Sr, Yb) can also act as 
ion sources due to photoionization or photoassociative 
ionization from their MOT beams I5ll. 
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FIG. 11. Plot of kinetic energy vs. time for Na + (orange) and 
Ca + (blue), (a) 10 Na + ions with only a Na background gas 
(Pback = 7 x 10" 9 torr), (b) 10 Na+ with 1 mK (10 -7 eV) 
MOT n = 5 x 10 10 cm -3 and r = 0.5 mm. (c) 10 Ca+ un- 
der the same MOT conditions. Ten ions can be cooled to a 
few Kelvin, but not to cold temperatures. 



The production of these extra NaJ and Na + ions pro- 
duced directly from the MOT results in an uncontrolled 
source of ions that can interfere with the controlled study 
of any other ions created within the hybrid trap (e.g., 
Na + initialized using REMPI). This extra co-trapped 
ion gas is much hotter than the MOT and has a strong 
(V oc 1/r) interaction with the Na + (or Ca + ) ions we are 
trying to cool. The thermalization of the different ion 
clouds would work against the ion- neutral (V oc 1/r 4 ) 
sympathetic cooling from the MOT. Furthermore, if ad- 
ditional ions are being created during the cooling period, 
the ion density may become too large and will eventually 
cause significant ion-ion rf heating. One can continuously 
and mass selectively quench the unwanted ions from the 
Paul trap via a well established experimental technique 
where an additional ac field is applied to either the end 
segments (on resonance with an ion's axial secular fre- 
quency) or the rf segments (on resonance with an ion's 
radial secular frequency) heating the ions above the pseu- 
dopotential trap depth. 

Unfortunately, these additional ac fields can have the 
side-effect of heating the ions we are trying to cool, de- 
spite the fact that the additional field is off-resonance 
with the cooled ion's mass dependent secular motion. 
Because of what we will call ac side-effect heating, we 
must use as low an ac field amplitude as possible. Addi- 
tionally, we found that radial quenching (as opposed to 
axial quenching) , as well as using higher harmonics of the 
quenched ion's secular motion helped reduce side-effect 
heating. This is likely because the radial trap depth is 
less than the axial trap depth and there is a larger differ- 
ence between various ion species' secular frequencies at 
higher harmonics. We found that the side-effect heating 
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FIG. 12. SIMION trajectories of two ions' equilibrating with 
the MOT (ion #1 in blue and ion #2 in black), (a) View 
along the axis of hybrid trap, (b) View from the side of the 
hybrid trap. The cooled ions initially result in a colder ion in 
the trap's center (ion #2) and a hotter ion (ion #1) in an orbit 
that is poorly overlapped with the MOT (rMOT = 0.5 mm). 



could therefore be significantly reduced (but not com- 
pletely removed) such that it offered effective quenching 
while negligibly increasing the equilibrium energy of the 
sympathetically cooled ions. 

While sympathetically cooling and trapping 10 Na + , 
we simulated the birth and simultaneous quenching of 
NaJ ions and found the process did not impede the cool- 
ing of the 10 Na + . Although an encouraging result, we 
should note that we could not simulate the actual Naj" 
birth rate (47j . because of computational limitations. 

We find a dramatic difference in the hybrid trap's abil- 
ity to sympathetically cool one ion (Fig. [3J compared to 
two or more ions (Fig. Hip . The main factor limiting the 
equilibrium energy of cooled multiple co-trapped ions is 
ion-ion rf heating, although atom-ion rf heating still ex- 
ists. In the presence of only a background gas, atom-ion 
heating contributes to a mean energy increase of 10 Na + 
seen in Fig. [TT] curve (a) [similar to single ion results 
in Fig. [3J curve (a)]. The heating is not due to ion-ion 
rf heating since (r rms ) is large enough that the ions are 
within the Mathieu regime [33[. 

Ten Na + (or Ca + ) ions [curve (b) and curve (c) of 
Fig. [TTJ respectively] cooled with a MOT density of 
n = 5 x 10 10 cm* 3 do not equilibrate at energies equiva- 
lent to sub- Kelvin temperatures (the ions are only cooled 
to energies equivalent to a few Kelvin due to ion-ion rf 
heating). However, the MOT sympathetic cooling should 
cause a significant extension in trapping lifetime since the 
ions are cooled well below the pseudopotential's radial 
trap depth of 0.94 eV and axial trap depth of 5 eV pjjj . 
Due to the ion-neutral mass ratio resulting in weaker 
atom-ion rf heating, Ca + equilibrates at a lower energy 
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FIG. 13. Plot of kinetic energy vs. time for 5 Ca 
ions. Curve (a) shows sympathetic cooling without crys- 
tallization (magenta) where Tmot = 100 fiK (10 -8 eV) and 
n = 5x 10 10 cm -3 . Curve (b) shows sympathetic cooling 
with crystallization (blue) n = 8 x 10 11 cm* 3 , where the im- 
age shows SIMIOn's rendering of 5 crystallized ions. 

than Na + . At this MOT density there was approximately 
no difference in the equilibrium energies for 2, 5, or 10 
sympathetically cooled ions. This is likely due the com- 
mon final energy barrier associated with the ion-ion rf 
heating. 

When trapping and cooling multiple ions a cold, nearly 
crystallized center was found with one or two hotter 
atoms orbiting around the periphery (see Fig. [T2"]) . At- 
tempts to improve overlap with the hotter orbiting ions 
by increasing the end segment voltage and placing a pos- 
itive bias on all four rf segments (effectively squeezing 
the cloud, i.e., increasing the overlap), did not signifi- 
cantly decrease the equilibrium energy. The lack of im- 
provement was consistent with the results discussed in 
the single ion case. 

2. Ion Crystal 

Decreasing the MOT temperature to 500 nK 
(6 x 10* 11 eV) slightly lowered the final energy of 
the ions but did not increase the cooling capacity enough 
to crystallize the entire ion cloud. Only a high density 
MOT (n > 1 x 10 11 cm* 3 ) can produce crystallization. 
Once cold enough to crystallize, we find a difference in 
final energies between the 2, 5, and 10 ion simulations. 
The minimum density needed for crystallization for 2 
ions with a 100 fiK (10* 8 eV) MOT was 4 x 10 11 cm* 3 , 
5 ions [shown in Fig. Q2] curve (b)] required at least 
8 x 10 11 cm* 3 , and 10 ions were never observed to 
crystallize, even at densities as high as 10 14 cm* 3 . 
When the MOT density is above the required minimum 
crystallization MOT density, further cooling can be 
realized. 
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FIG. 14. Plot of (r rms ) vs. 5 ions average kinetic energy for 
two MOT densities n = 5 x 10 10 cm" 3 
which does not crystallize and n ; 



(green or light gray) 



8 x 10 10 



(black) 



which does crystallize. Marker (1) denotes the ion-ion rf heat- 
ing barrier, that is not breached at the lower MOT density. 
At higher MOT density crystallization is reached at marker 
(2). 



One can determine the equilibrium ion crystal spacing 

25 /urn 



(Mpa 3 ) 



1/3 



(11) 



of two ions by equating the restoring force due to the 
trap 's axial potential and the ion-ion Coulomb repulsion 
33]. In Eq. (fTTj) M is the atomic mass of one ion in 
atomic mass units and / is the rf driving frequency in 
MHz. The simulated 2 and 5 ion equilibrium separation 
of the crystal shown in Fig. [13] agrees within a few mi- 
crometers with the value obtained using Eq. (jlip . 

To support the claim that ion-ion rf heating is the 
mechanism that determines the final energy for sympa- 
thetic cooling of multiple ions, we examined the correla- 
tion between mean energy of the cooled ions and the ion 
cloud's (r rms ) (see Fig. [14]). While at initially large (r rms ) 
and kinetic energy (i.e., within the Mathieu regime) there 
is little difference in the cooling by the high or low den- 
sity MOT, except for small fluctuation due to atom-ion 
rf heating. As the (r rms ) decreases, we begin to enter the 
chaotic regime; the heating rate begins to fluctuate as a 
function of the ion cloud (r rms ) resulting in both a loss of 
clear position and energy correlation [shown by marker 
(1) in Fig. [14]. This could be thought of as an energy 
~ (?*rms) barrier. Only with the cooling capacity of the 
higher MOT density can the ion cloud move past the ion- 
ion rf heating barrier into an ion-crystal phase. Once in 
the approximately constant (r rms ) crystal state [shown 
by marker (2) in Fig. Q3] the ions can then be cooled 
further (by reducing the small oscillation amplitudes). 

To test the effect of atom-ion rf heating on crystal- 



lization, we simulated the sympathetic cooling of 5 ions 
that were more massive than Ca + (mj/m n ~ 7.52), but 
assumed the same elastic scattering rate coefficient as 
that of Na-Ca + . We found that these more massive ions 
cool to a lower final energy than Ca + (as can be expected 
with reduced atom- ion rf heating) , but that the minimum 
MOT density required to crystallize the ions is approxi- 
mately the same. Hence, the only way to achieve cold or 
ultracold ion cloud temperatures is to have a high enough 
MOT density to overcome the ion-ion rf heating. 



IV. CONCLUSION 

We simulated sympathetic cooling of a single ion and 
multiple ions (2 < Nj on < 10) in a hybrid trap. Our 
findings demonstrate that a MOT with a low density 
~ 10 9 cm -3 and modest 1 mK MOT temperature can 
cool a single ion to ultracold energies within seconds, even 
in instances of equal ion and neutral mass. Therefore, a 
BEC is not required to achieve sympathetic cooling of a 
single ion in a hybrid trap. 

To achieve the most effective cooling, we found that 
it is critical that the MOT be concentric with the ion 
cloud and as dense as possible. The MOT cooling rate is 
larger than the atom-ion rf heating rate for only part of 
the full range of stable qt values. Decreasing the MOT 
temperature does decrease the final ion energy. However, 
it does not do so indefinitely, since the atom-ion rf heating 
rate also increases. 

Modest MOT conditions can also sympathetically cool 
more than one trapped ion, although not to sub-Kelvin 
temperatures. High MOT densities (n > 1 x 10 11 cm~ 3 ) 
or BEC densities are needed to overcome the ion-ion rf 
heating, crystallize the ions, and allow for the possibility 
of further cooling toward ultracold temperatures. How- 
ever, this appears to only be experimentally feasible for 
a small number of ions (Nj on < 10). 

The initial ion cloud's temperature can be determined 
via simulation, as it depends primarily on the initial root- 
mean-squared position of the ion cloud (r rms ) rather than 
the number of ions created or the temperature of the 
neutral gas from which they are born. Last, it should be 
possible to perform mass selective ion quenching of one 
species without significantly heating other ions that are 
being sympathetically cooled. 
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FIG. 15. ID Atom-Ion rf heating phase space diagrams [2ll |. 
(a) ion phase space trajectory following simple harmonic mo- 
tion at instance of rf phase = 0. (b) ion phase space tra- 
jectory at instance of rf phase 7^ 0. Arrows 1-3 and l'-4' 
(beginning with a solid dot and ending with an arrow head) 
show various instantaneous speed changes caused by atom-ion 
collisions, i.e., the energy gained or lost in a collision results 
in a new path in phase space. The thick (black) ellipse is the 
pre-collision path, while a hotter end path corresponds to a 
larger ellipse (red) and a colder one is a smaller ellipse (blue or 
magenta). Atom-ion rf heating is shown in arrow 3', where an 
instantaneous speed decrease results in a larger phase space 
ellipse, and therefore an increase in total energy. 



Appendix A: Paul Trap Heating Mechanisms 

Linear Paul traps are susceptible to a variety of pos- 
sible inherent heating mechanisms: instantaneous col- 
lisional heating, two types of rf heating (atom-ion and 
ion- ion heating), excess micromotion heating, and insta- 
bility heating. We wish to emphasize that each of the 
five mechanisms mentioned are unique and it is best not 
to group them all together under the imprecise heading 
of rf heating. 

Both instantaneous collisional heating and atom-ion rf 
heating within a buffer gas have been clearly explained 
by S. Schwarz using phase space diagrams 21] ; a modi- 
fied version of his explanation is shown in Fig. [TOJ Each 
phase space diagram shows the one dimensional (ID) 
momentum-position relationship of a single ion's secular 
motion in a Paul trap during a specific instance of the rf 
driving field's phase (i.e., at a single instance during the 
micromotion). The total energy of the ion (kinetic plus 
potential) is equal to the area of the phase space ellipse. 
At any instant in time, the ion has a distinct position 



in phase space constrained to a point on a closed path. 
Simultaneously, as the phase of the rf potential evolves in 
time, the closed phase space path rotates and stretches. 

We will assume that the interaction time during a col- 
lision between an ion and atom is much shorter than 
the period of the ion's secular motion or micromotion. 
Therefore, collisions are represented as vertical arrows 
in Fig. 1151 where each collision event (beginning with a 
solid dot and ending with an arrow head) instantaneously 
moves the ion to a new path in phase space consistent 
with the ion's new post-collision total energy. 

A collision causing an instantaneous speed increase, 
represented by arrows 1 and 1', always results in a total 
energy increase (independent of rf phase). We will de- 
fine this event as instantaneous collisional heating. For 
example, this type of collision might occur because the 
ion had a head-on collision with a hotter background gas 
atom, because it was hit by an atom traveling in the same 
direction with any non-zero amount of energy, or because 
it was hit by an atom with a mass ratio mj /m n < 1 [l8j . 
Instead, if the ion undergoes an instantaneous speed de- 
crease, depicted by all other arrows in Fig. [TOJ the ion 
may or may not result in a total energy decrease (de- 
pending on the instance of the rf phase at the moment 
of the collision) . 

When the ion follows a simple ID harmonic phase 
space ellipse at the rf phase = instant, any instanta- 
neous speed decrease will always result in a lower en- 
ergy phase space ellipse and therefore instantaneous col- 
lisional cooling. These collisions are illustrated by arrows 
2 and 3 in Fig.fTBTa). which show a small and large speed 
change, respectively. The small speed change collisions 
arise when mi/m n > 1, while the large ones occur when 
mi/m n w 1. 

Atom-ion rf heating results from collisions that oc- 
cur during an instance when the rf phase does not fol- 
low a simple harmonic phase space path [phase =/= 
in Fig, froth)]. Collisions causing a large instantaneous 
speed decrease (dramatically interrupting the ion's mi- 
cromotion) can result in a total energy increase [see ar- 
row 3'], i.e. atom- ion rf heating. For the case where 
mi/m n > 1, a head-on collision may result in a smaller 
instantaneous speed decrease and less atom- ion rf heating 
[non- heating collisions depicted by arrow 2'] [l8|, [2l[ . 

The MOT in our hybrid trap has a high density over 
a finite region at the center of the Paul trap. Therefore, 
all collisions with the MOT, including head on collisions 
with mi /m n « 1, result in little to no atom-ion rf heating 
until the ion's secular motion amplitude is smaller than 
the MOT radius tmot (see arrow 4') where the micro- 
motion amplitude is smaller 10]. For buffer gas cooling, 
atom-ion collisions can occur throughout the entire vol- 
ume of the Paul trap; hence cooling ions with neutrals of 
equal mass results in no net temperature change. 

Our discussion thus far has assumed head-on collisions 
in ID. In 3D, the majority of collisions will be glancing, 
resulting in smaller changes in ion speed than predicted 
by the ID model. Also, the colliding MOT atom's veloc- 
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FIG. 16. Sketch of ion-ion rf heating (green), heating rate 
vs. (r rms ) At large (r rms ) an ion cloud is unaffected by rf 
heating and the ions undergo motion described by the Math- 
ieu equation. As the ions are cooled and {r rms ) decreases, the 
Coulomb interactions cause chaotic motion among the ions 
and absorption of the rf driving field's energy, i.e., ion-ion rf 
heating. If cooled further, the ion-ion rf heating rate eventu- 
ally drops with decreasing (r rms ) and the ion cloud becomes 
an ion crystal. 



ity in the laboratory frame is small compared to the ion's 
velocity during portions of the ion's secular motion, but 
is not zero velocity at all times. Therefore, most equal 
mass head-on collisions do not force the ion to come to 
a dead stop (depicted by arrow 3'). This suggests that 
even when the ion's amplitude is smaller than the MOT 
radius further cooling is possible for equally massive ion 
and neutral species. 

The heating mechanism known as ion-ion rf heating 
was first explained by R. Bliimel et al. and we present a 
sketch of their main results in Fig. [16] to briefly illustrate 
this effect P, EH . 

For simplification, let us temporarily assume there is 
no atom-ion rf heating. The root-mean-squared posi- 
tion of the ion cloud, denoted by (r rms ) is given by 



i_ 



where r; is the position of the Ith ion 



in an ion cloud containing N ions. The ions are initially 
far enough apart that they do not experience any heat- 
ing. As these ions are cooled (by any means) the axial 



amplitude decreases [see Eq. therefore the (r rms ) ini- 
tially decreases without changing the heating rate. Since 
the motion of individual ion's can still be described by 
the Mathieu equation, this is what Bliimel called the 
"Mathieu regime" and can be seen in Fig. [TBI Ion-ion 
rf heating begins once the (r rms ) becomes small enough 
that Coulomb interactions between ions perturb the ions' 
motion, resulting in deterministic chaos [52j. The chaos 
smears the well defined frequency spectrum indicative of 
the Mathieu regime, allowing for resonant absorption of 
the rf driving field's energy which causes the ions to heat 
up. Bliimel called this region the "chaotic regime." If 
the cooling capacity is great enough to overcome the rf 
heating then a phase transition can occur, resulting in an 
ion crystal (the ions lie in a chain along the trap axis). 
At the boundary between crystallization and the chaotic 
regime, ions undergo quasiperiodic motion which Bliimel 
called the "quasiperiodic regime." Only after getting be- 
yond the ion- ion rf heating hurdle can crystallization and 
continued ion cooling occur. 

Excess micromotion occurs when an ion is displaced 
away from the quadrupole field's nodal line (found along 
the X3 axis of the trap). Experimentally, the displace- 
ment can occur inadvertently from stray electric fields 
within the laboratory, or be actively controlled by ei- 
ther adding dc potentials (dipole configuration) to the rf 
segments or by changing the phase relationship between 
each diagonal pair of rf segments away from the ideal 
180° [1, yij]. Once displaced, the ion will have an exces- 
sive or larger micromotion amplitude compared to when 
it is not displaced, thereby increasing the ion's mean en- 
ergy. Since the mechanism is driven by the rf field, sym- 
pathetic cooling cannot reduce excess micromotion [34j . 
This mechanism is different from rf heating, as it can oc- 
cur in the absence of any ion-ion or atom-ion collisions. 

Instability heating is strictly a multiple-ion phe- 
nomenon, caused by the random interaction between ions 
resulting in instantaneous unstable qi values. Unlike ion- 
ion rf heating, this type of heating can be eliminated by 
choosing a low qi value, which can also reduce rf heating 
(but cannot eliminate it) (l7| . As we do not simulate 
nodal line displacement fields and we use low values, 
both excess micromotion heating and instability heating 
can be ruled out as having an effect within the context 
of the simulations presented in this paper. 
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